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Abstract 

The Gribov ambiguity problem is studied for compact U{1) lattice theory within 
the Lorentz gauge. In the Coulomb phase, it is shown that apart from double Dirac 
sheets all gauge (i.e. Gribov) copies originate mainly from the zero-momentum 
modes of the gauge fields. The removal of the zero-momentum modes turns out to 
be necessary for reaching the absolute maximum of the gauge functional F{6) . A 
new gauge fixing procedure - zero-momentum Lorentz gauge - is proposed. 

1 Introduction 

To gain a better understanding of the structure of the lattice theory and to interprete 
correctly the numbers obtained in Monte Carlo simulations, it is very instructive to 
compare gauge variant quantities such as gauge and fermion field propagators with 
the corresponding analytical perturbative results. In this respect, compact U{1) 
pure gauge theory within the Coulomb phase serves as a very useful 'test ground', 
because in the weak coupling limit this theory is supposed to describe noninteracting 
photons. However, previous lattice studies |jl|, |^, ||, ^, |^ have revealed some rather 
nontrivial effects. It has been shown that the standard Lorentz (or Landau) gauge 
fixing procedure leads to a r-dependence of the transverse gauge field correlator 
^t{t]P) being inconsistent with the expected zero- mass behavior [^. Numerical 
0, H, H and analytical studies have shown that there is a connection between 
'bad' gauge (or Gribov) copies and the appearance of periodically closed double 
Dirac sheets (DDS). The removal of DDS restores the correct perturbative behavior 
of the transverse photon correlator r2-(r;p) with momentum . However, it 



does not resolve the Gribov ambiguity problem completely. Other Gribov copies 
connected with zero-momentum modes of the gauge fields still appear, which can 
'damage' such observables as the zero-momentum gauge field correlator r(r; 0) or 
the fermion propagator r^(r) ||, 0. 

After having understood, why gauge variant lattice correlators behave unex- 
pectedly from the perturbation theory point of view, one can search for the 'true' 
gauge fixing procedure. This constitutes the main goal of this note. We propose 
a zero-momentum Lorentz gauge (ZML), which permits to get rid of the lattice 
artifacts and provides correct values for various correlation functions. 

We are going to compare also the standard Lorentz gauge fixing procedure (LG) 
P with the axial Lorentz gauge (ALG) proposed in [|^. We show that ALG 
produces just the same problems as LG does and, therefore, cannot resolve the 
Gribov ambiguity problem. 



2 Gauge fixing : Lorentz gauge and axial Lorentz 
gauge 

The standard Wilson action with f/(l) gauge group is [|ll| 

-^If/) =/5EE(l- cos ^.,^.) , (2.1) 

where the link variables are = exp^iO^fj) G U{1) and 9^^ G (— 7r,7r] . The 
plaquette angles are given by 6^.^,^ = 6^.^ + 6x+ii-u - O^+t,-^ - 6^^ . This action 
makes the part of the full QED action Sqed, which is supposed to be compact if 
we consider QED as arising from a subgroup of a non-abelian (e.g., grand unified) 
gauge theory fl^ . 

The plaquette angle 6p = 6x-fj,u can be split up: 6p = [6p] + 2iTnp, where 
[6p] G (— tt; tt] and np = 0,±l,ib2. The plaquettes with np are called Dirac 
plaquettes. The dual integer valued plaquettes m^^^i, = \e^ypanx^pcj form Dirac 
sheets [|13| . In our calculations we monitored the total number of Dirac plaquettes 



N^^^ for every plane (/x; v) = (1; 1), ... , (3; 4) and N^p = max(^,^) N^^^p . The 
appearance of periodically closed double Dirac sheet means that at least for one 
out of six planes (/xz/) the number of Dirac plaquettes N^^'^ should be 

> 2^ . (2.2) 

For example, on the lattice 12 x 6^ the appearance of DDS means Npip > 72. 
In lattice calculations the usual choice of the Lorentz (or Landau) gauge is 

4 

J2d^sm9xp = 0, (2.3) 
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which is equivalent to finding an extremum of the functional F[6) 

= T7E^-W ; = iEfcos^,^ + cose,._^;^ 

with respect to the (local) gauge transformations 



(2.4) 



a 



Xfl 



A U A* - 

^^x'-^Xfl^^X+fl 



A^ = exp{in^} e f/(l) . 



(2.5) 



The standard gauge fixing procedure (referred in what follows as LG) consists 
of the maximization of the value Fx{6) at some site x under the local gauge 
transformations A^^ , then at another site and so on. After a certain number of 
gauge fixing sweeps a local maximum Fmax of the functional F{9^) is reached. In 
order to improve the convergence of the iterative gauge fixing procedure one can use 



the optimized overrelaxation procedure |]I4[ with some parameter a which depends 
on the volume and (3 . For convergence criteria we use 



max 



E sin 0, 



XfJ, 
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V ^ 



E sin 6*. 



< 10" 



In Figure |l]a we show a typical Monte Carlo time history of F^ax obtained using 
the standard Lorentz gauge fixing procedure. Partially, the comparatively big dis- 
persion of Fmax is due to periodically closed double Dirac sheets 0, |^ . In Figure 
|T|b we show the corresponding time history of N^p = max(^,^) Nj!^p . For the given 
lattice size and /3-value ~ 20% of configurations turn out to possess DDS. Exactly 
because of the configurations containing DDS the transverse photon correlator 



rr(r;p) = (<l>(r;p) $*(0;p)) , <I>(r;p) = Y.exp{ipx + '-p,,) sin^,^,^ (2.6) 

X 

( = 1,3, p = (0,p, 0) ) exhibits an unphysical 'tachion-like' behavior p|, 0, ||. 
In Figure ^ we show the normalized photon correlator Tt{t] p) /rT(0; p) for lowest 
non- vanishing momentum and LG together with results obtained with other Lorentz 
gauge fixing procedures to be discussed lateron. We clearly see for LG the deviation 
from the expected zero-mass behavior. All the observations described above have 
not been seen changing, when j3 and/or the lattice size were increased considerably 

i- 

It is a long-standing believe [|l5l that the 'true' gauge copy corresponds to the 



absolute maximum of F{6). Thus it would be highly desirable to find a Lorentz 
gauge fixing procedure which produces a (unique) gauge copy of the given configu- 
ration with the absolute maximum of F{9) . However, if one repeatedly subjects 
a configuration {Ox^i} to a random gauge transformation as in Eq. (|2.5|) and then 
subsequently applies to it the LG procedure, one usually obtains gauge (Gribov) 
copies with different values of F„ 



max- 
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An attempt to resolve this problem has been made in where a modified 
Lorentz gauge fixing procedure has been proposed. This procedure (which we refer 
to as the axial Lorentz gauge fixing or ALG) consists of the two steps: 



i) first transform every configuration to satisfy a maximal tree temporal gauge 
condition ('axial' gauge) |T6| ; 

ii) then apply the Lorentz gauge fixing procedure.^ 

An axial gauge with a chosen maximal tree is unique by definition. In practice, 
this is easily checked by random gauge transformations applied first. Consecutive 
gauge fixing steps - e.g. the Lorentz gauge iterations - will lead always to the 
same result as long as we do not change the detailed prescription for these steps. 
The question is, whether this 'unique' Lorentz gauge obtainable for each gauge field 
configuration resolves the problems mentioned above. The answer is 'no'. We do 
not find the absolute maximum of the gauge functional in the majority of the cases. 
There is a quite high percentage of Gribov copies left containing DDS (around 
10% for (3 = 1.1 and a 12 x 6'^ lattice). As a consequence the transverse 
non-zero momentum photon propagator does not come out correctly again. The 
corresponding data are shown in Figure 0, too. Doubling of the linear lattice size 
and enlarging /3 (we checked j3 = 2.0 ) do not improve the behavior. 

In paper p| a Lorentz gauge fixing prescription with a preconditioning step 
based on a non-periodic gauge transformation was proposed. The latter has been 
chosen in such a way that the spatial Polyakov loop averages were transformed into 
real numbers as a first step. We convinced ourselves that the Lorentz gauge fixed 
configurations with very high probability did not contain DDS. As a consequence, 
the photon correlator becomes correct. However, this gauge fixing procedure does 
not provide the absolute maximum of the Lorentz gauge functional, too. Thus, the 
lattice Gribov problem for QED in the Coulomb phase has not been solved. 



3 Zero— momentum Lorentz gauge (ZML) 

As is well known, apart from the local symmetry Wilson action S has an additional 
(global) symmetry with respect to non-periodic transformations 

or equivalently 



^In an analytically known local solution of the Lorentz gauge condition has been applied, 
instead of an iteratively obtained one. But, as we have convinced ourselves, this solution does not 
resolve the problems we are discussing here. 
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Note that these transformations do not spoil the periodicity of the gauge fields Ux^ 
and 6^^ , respectively. The transformation in Eq.( p.2|) changes the zero-momentum 
mode (p^ of the link angle 9^^ = (p^ + dO^^ , where ^^x^t = , and a proper 
choice of the parameters can make 0^ equal to zero. 

It is rather evident that for the infinitesimal fluctuations 5d.j;^ about the zero- 
momentum mode 0^ the absolute maximum of the functional F{9) corresponds 
to the case 0^ = . 

Our main statement (to be proved in this section) is that in most of the cases 
(^99.99%) gauge copies of the given configuration O^fj. are due to a) double Dirac 
sheets; h) the zero-momentum modes of this field. It is the exclusion of DDS and 
of the zero-momentum modes that permits to obtain a gauge copy of the given 
configuration with the absolute maximum of the functional F{6) . 

We define our gauge fixing prescription (which we refer to as the zero-momentum 
Lorentz gauge or ZML gauge) as follows. Every iteration consists of one sweep with 
(global) transformations as in Eq. ( p.2|) and one sweep with local gauge transforma- 
tions as in the standard Lorentz gauge fixing procedure. 

In Figures |^,b we show time histories of Fmax and of Nnp for the gauge fixing 
procedure with the minimization of the zero-momentum modes of the gauge field. 
Both the time histories should be compared with the corresponding ones for the 
standard Lorentz gauge shown in Figures |l]a,b. One can see that after suppression 
of the zero-momentum modes the average value of Fmax became essentially larger 
as compared with the LG case. At the same time the number of DDS (to be 
precise, the number of configurations with Nop > 72) has drastically decreased 
(~ 1%). Typically, those configurations with DDS have smaller values of F^ax 
than configurations without DDS have. 

It is a very easy task to remove the remaining DDS. If a DDS in a ZML gauge 
has really appeared, then perform a random gauge transformation to the same 
field configuration and repeat the ZML procedure again. As a result gauge field 
configurations do not contain DDS. 

To convince ourselves that the ZML gauge fixing procedure provides an absolute 
maximum of the functional F{9) we generated many random gauge copies for every 
thermalized configuration {Ox^i} ■ The number of these random gauge copies N^c 
varied between 10 and 1000 for different /5's and lattices. Let {^'i:'^} 
the j*^ gauge copy of the configuration obtained with the random gauge 

transformation as in Eq.(p.5|), j = 1, . . . , Njic ■ For any configuration \^0xfi\ "we 
define a 'variance' SFmax{0) of Fmax{d) 

6Fmax{0) = max(F„,,(e«) - Fmax{0^'^)) ; i, j = I, . . . , Nrc ■ (3.3) 

Different gauge copies can, in principle, have different values of F^axi^) ■ Therefore, 
its 'variance' SFmaxi^) can be non-zero. For example, the standard gauge fixing 
procedure (LG) described above gives typical values SF^axi^^) ~ 0.06 0.07 for 
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Nrc = 10 . As an example we show in Figure |a a time history of the 'variance' 
6Ffnax{(^) ioT the standard LG and for ZML. In the case of ZML this 'variance' is 
zero 0. 

In Figure ^3 we compare values of F^ax for these two gauges. The lower bro- 
ken line for the LG corresponds to the first gauge copy which is just a thermalized 
configuration produced by our updating subroutine. Applying random gauge trans- 
formations we produced different random gauge copies {Nuc = 10 in this case) and 
then chose the gauge copy with the maximal value of F^ax ■ The corresponding 
values are represented by the solid LG line in Figure |^. One can see that even 
this is far below the solid line which corresponds to the ZML gauge. Increasing the 
number Njic up to 1000 does not change this conclusion. 

One reason for the non-zero value of the 'variance' is the appearance of gauge 
copies containing double Dirac sheets 0, H, |^. Gauge copies with double Dirac 
sheets have typically much lower values of F^axi^) as compared to Fmax for gauge 
copies without DDS. Moreover, a contribution of configurations with DDS 'spoils' 
the photon correlator rT{T;p) with p ^ which leads to a wrong dispersion 
relation inconsistent with the dispersion relation for the massless photon. Double 
Dirac sheets represent a spectacular example of lattice artifacts which can lead to 
a misleading interpretation of the results of numerical calculations (see also ||17|| ). 

But, as we have seen, the mere exclusion of gauge copies with DDS in the 
standard Lorentz gauge fixing procedure does not yet provide a zero value of the 
'variance' SFmaxi^) ■ Different gauge copies (without DDS) still have different 
values of Fmax{d) ■ The exclusion of the zero- momentum modes turns out to be 
sufficient to remove the ambiguity. 

4 Conclusions 

Now let us summarize our findings. 

Our main result is that ~ 99.99% (if not all) gauge copies for the given gauge 
field configuration are due to two reasons : 

a) periodically closed double Dirac sheets; 

2) the zero-momentum mode of the gauge field 9^^ ■ 

We didn't find any other reason for the appearance of the gauge copies. 
The minimization of the zero-momentum mode can be performed sweep by sweep 
using a global transformation as in Eq.( p.2| ). We proposed a modified gauge fixing 
procedure consisting of local gauge transformations in Eg. ( |2.5|) and global trans- 
formations in Eq.( p.2| ) (ZML gauge). The exclusion of the DDS (which appear 

In rather rare cases (<0.01%) in ZML gauge 5F^ax{0) 7^ . So far we have no explanation 
for it. 
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rather rarely in the ZML gauge) can be easily performed on the algorithmic level 
as described in the text. 

The application of the ZML gauge fixing procedure provides us with absolute 
maximum of the functional F{6) . 

We have shown that the gauge fixing procedure with axial gauge preconditioning 
(ALG) cannot solve the problem of the Gribov ambiguity in this theory. The axial 
gauge preconditioning cannot exclude the appearance of DDS as well as of nonzero 
values of the zero-momentum mode. 

In this paper we present our results only for the bosonic sector of the theory. 
However, we believe that this study solves ultimatively the Gribov ambiguity puzzle 
in the case of quenched compact QED within the gauge as well as the fermionic 
sector. Work on the fermion case is in progress 0. 

The inclusion of dynamical fermions changes somewhat the symmetry group of 
the full action, i.e. the parameters in Eq.'s (|3.1|), ( p.2[ ) can have only discrete 
values. This case needs some additional study. 

Acknowledgements 

Financial support from grant INTAS-96-370, RFBR grant 99-01-01230 and the 
JINR Heisenberg-Landau program is kindly acknowledged. 

References 

[1] A. Nakamura and M. Plewnia, Phys. Lett. 255B (1991) 274. 

[2] V.G. Bornyakov, V.K. Mitrjushkin, M. Miiller-Preussker and F. Pahl, 
Phys. Lett. 317B (1993) 596. 

[3] V.K. Mitrjushkin, Phys. Lett. 389B (1996) 713. 

[4] V.K. Mitrjushkin, Phys. Lett. 390B (1997) 293. 

[5] l.L. Bogolubsky, L. Del-Debbio and V.K. Mitrjushkin, JINR E2-99-43, 
SHEP/99/01, |hep-lat/99030T5| . 

[6] Ph. de Forcrand and J.E. Hetrick, Nucl. Phys. B (Proc.Suppl.) 42 (1995) 861. 

[7] l.L. Bogolubsky, A. Hoferichter, V.K. Mitrjushkin, M. Miiller-Preussker and 
P. Peter, in progress. 

[8] V.N. Gribov, Nucl. Phys. B139 (1978) 1. 

[9] J.E. Mandula and M. Ogilvie, Phys. Lett. 185B (1987) 127. 



7 



[10] A. Nakamura and R. Sinclair, Phys. Lett. 243B (1990) 396. 
[11] K. Wilson, Phys. Rev. DIO (1974) 2445. 

[12] A. M. Polyakov, Gauge fields and strings, Harwood Academic Publishers, 
(1987). 

[13] T.A. DeGrand and D. Toussaint, Phys. Rev. D22 (1980) 2478. 
[14] J.E. Mandula and M. Ogilvie, Phys. Lett. 248B (1990) 156. 

[15] D. Zwanziger, Nucl. Phys. B364 (1991) 127; B378 (1992) 525. 
[16] M. Creutz, Phys. Rev. D15 (1977) 1128. 

[17] V. Grosch, K. Jansen, T. Jersak, C.B. Lang, T. Neuhaus and C. Rebbi, 
Phys. Lett. 162B (1985) 171. 



Figure captions 

Figure 1. Time history of F^ax (a) and Ndp (b) at /3 = 1.1 on the 12 • 6^ lattice 
in the standard Lorentz gauge. 

Figure 2. Transverse propagator at /3 = 1.1 on the 12-6^ lattice in three different 
gauges. 

Figure 3. Time history of -F^ax (a) and Ndp (b) at /3 = 1.1 on the 12 • 6^ lattice 
in ZML gauge. 

Figure 4. Time history of the 'variance' SF^^ax (a) and F^^ax (b) at /3 = 1.1 on 
the 12 • 6^ lattice. LG solid line corresponds to N^c — 10. Gauge copies with DDS 
are excluded. 
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Figure 1: Time history of Fmax (a) and Ndp (b) at /9 = 1.1 on the 12 • 6^ lattice 
in the standard Lorentz gauge. 
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Figure 3: Time history of Fmax (a) and Ndp (b) at /9 = 1.1 on the 12 • 6^ lattice 
in ZML gauge. 
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Figure 4: Time history of the 'variance' 5Fmax (a) and Fmax (b) at /3 = 1.1 on the 
12 • 6^ lattice. LG sohd hne corresponds to Nrc — 10. Gauge copies with DDS are 
excluded. 
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